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INTRODUCTION 


The ion beam divergence characteristics of ion accelerator systems 
have been -studied by many workers in recent years This increased 
interest in ion-optical phenomena is a result of the direct application 
of ion beams to three new technology areas that have evolved rapidly in 
the last decade. Briefly, these areas are the following. Electric 
propulsion of space vehicles; where ion thrusters provide thrust by 
ejecting beams of energetic ions. Neutral beam heating of fusion 
plasmas; where deuterium ions are accelerated, resonant charge ex- 
changed, and injected into a fusion device as neutral particles to heat 
directly the confined plasma. Ion beam sputtering, deposition and 
milling; where carefully controlled ion fluxes are used to alter the 
surface structure and/or composition of thin films for semi-conductor 
and other applications. Well focused, high current, ion beams are a 
necessary requirement in each of the aforementioned areas. 

Presently, the ion acceleration and focusing process is fairly 
well understood. Using the data available in Refs. 6 and 7, one can 
quite accurately predict the ion beam divergence characteristics of an 
arbitrary two or three-grid ion accelerator system design over a wide 
range of specified operating conditions. However, the process by which 
ions are extracted from the discharge plasma, prior to their accelera- 
tion and focusing by the accelerator system, is not well understood. 

It is known that a plasma sheath is formed and that this sheath is a 
transition region separating the discharge plasma (i.e., the ion pro- 
duction region where local ion and electron number densities are equal) 
from the ion acceleration region where ions only are present. This 
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plasma sheath is expected to vary in position and shape as e result of 
plasma density (i.e., beam current), accelerator system geometry and 
accelerator system potential variations. Knowledge of how these para- 
meters affect the physical characteristics of the plasma sheath would 
advance the understanding of the role this sheath plays in ion optical 
phenomena. These physical characteristics include, the divergence of 
the beam ions ejected from the accelerator system, the efficiency with 
which ions are extracted from the discharge plasma, the ion bombardment 
and erosion of various accelerator system components and the limit to 
the ion current which can be extracted from the plasma. 

. Some theoretical and experimental work has been done to determine 

the physical characteristics of the plasma sheath transition region 

2 4 5 8—1 0 

separating the discharge plasma and ion accelerator system. » » » “ 

However, the theoretical model predictions are somewhat conjectural 

since they rely on various initial assumptions which may not necessarily 

be valid in the actual physical situation. This suspicion is rein- 

3 6 7 

forced with earlier work by the author * ’ showing serious discrep- 
ancies between theoretically predicted and experimentally observed ion 
beam divergence characteristics of various ion accelerator systems. 

I • 

These discrepancies are believed to result from the inability of the 
theoretical models to describe the plasma sheath transition region 
adequately. Similarly, the available experimental data is limited to 
a couple of idealized cases; which are photographic in nature and 
present the plasma sheath as a discontinuity and are of limited use. 
Therefore, a thorough experimental investigation of the physical 
structure and behavior of the plasma sheath adjacent to an ion accel- 
erator system was undertaken. The basic physical phenomena uncovered 
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would provide a substantially better understanding of the ion extrac 
tion and pre-ion acceleration process from a plasma and aid in the 
development of improved theoretical models. 


ACCELERATOR SYSTEM SCALING 


Ion Beam Formation 

Most ion accelerator systems consist of either two or three grids 

(i.e., electrodes) with each grid perforated by numerous holes, usually 

in a hexagonal array. These grids are electrically isolated from each 

other and orientated so that their holes are aligned. The assembled 

accelerator system is positioned adjacent to a discharge chamber capabl 

* 

of producing the desired plasma conditions. Ion acceleration is 
achieved by applying a potential distribution to the accelerator system 
such that plasma electrons are repulsed and plasma ions are accelerated 
by the applied electrostatic body forces. Figure la depicts a portion 
of a multi -aperture two-grid ion accelerator system showing the coaxial 
hole geometry and ion beamlet formation. Figure lb illustrates the 
variation in electrical potential associated with the grid geometry of 
Fig. la. As shown in Fig. lb, the discharge chamber is at a potential 
a few tens of volts positive of the screen grid which is at a high 
positive potential. The purpose of the screen grid is to prevent di- 
rect access of the accelerated ions to the accelerator grid webbing. 

It is the accelerator grid, with its high negative potential, that 
actually accelerates the ions from the discharge plasma. A third, 
decelerator grid, may be positioned after the accelerator grid. This 

grid is usually maintained at facility ground potential and can reduce 

4 7 

ion beam divergence significantly. ’ 


* 


Here, as in the following sections, positive ions only are being 
considered. 


POTENTIAL 
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SCREEN HOLE SCREEN 
PLASMA SHEATHv GRID 


ACCELERATOR 

GRID 



Fig. 1 a) Accelerator system geometry. 

b) Accelerator system potential variation. 
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The simple two-grid ion accelerator system shown in Fig. 1 
embodies those principles of operation of most importance to the physi- 
cal phenomena occurring in the screen hole plasma sheath transition 
region. It is in this region that the plasma electrons begin to feel 
the effect of the negative accelerator grid and are retarded. Since 
the plasma electrons have a distribution of energies, their penetration 
depth into the ion accelerator system is varied; giving rise to the 
large extent of screen hole plasma sheath depicted in Fig. la. 

Geometric Similarity 

In order to make accurate measurements of the physical character- 
istrics of the screen hole plasma sheath (i.e., position, shape and 
structure), it was necessary to use a large accelerator system so ade- 
quate spatial resolution could be achieved. This was done by increasing 
the screen hole diameter from a value of 2.0 mm, which is typical of 
ion thruster and ground application ion source accelerator system de- 
signs and where the screen hole plasma sheath width was estimated to 
be 0.3 mm, to a diameter of 12.7 mm where a plasma sheath width of 
'v 2.0 mm was estimated. The other geometrical accelerator system 
parameters were also scaled directly by the ratio of these two diam- 
eters so geometric similarity was maintained. It was important to 
determine whether the maximum normalized perveance per hole^ and ion 
beam divergence characteristics of such large accelerator systens were 
any different from the much smaller accelerator systems typical of most 
ion sources. Careful probing of the ion beam emerging from the large 
accelerator systems showed the divergence characteristics of these grid 
sets were identical, within experimental error ±1.0°), to those of 
the smaller accelerator systems. Typical results are presented in 
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Table 1 where the ion beam divergence angle has been defined as the 
cone half angle, a, enclosing 955S of the total 'beam current.^ 

TABLE 1. Large Accelerator System Ion 

Beam Divergence Characteristics. 


NP/H , 
(amp/volt3/2) 

a, measured 
(degrees) 

a, predicted 
(degrees) 

1.08x10-^ 

16.0 

15.9 

2.17x10-9 

11 .3 

12.4 

3.25x10-9 

11 .5 

12.5 


Figure 2 shows the effect of screen hole diameter variations on 
the maximum normalized perveance per hole for grid sets having similar 
nondimens ional i zed geometric parameters. Here, the total accelerating 
voltage Vy, the discharge voltage and the net-to-total accelerating 
voltage ratio R were held constant at the values indicated. These 
parameters were not dependent on the accelerator system geometry scale 
changes. From Fig. 2 it is evident that there is little effect on the 
maximum normalized perveance per hole as the screen hole diameter is 
increased above 2.0 mm. The significance of this result can be fully 
appreciated when it is understood that a screen hole diameter change 
of 2.0 to 12.7 mm corresponded to a factor of fourty decrease in beam 
current density. It is not important here, but some effects have been 
found for smaller screen hole diameters than 2.0 mm.® 

In summary, scaling up ion accelerator systems to larger sizes 
does not appear to affect ion beam divergence characteristics or the 
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maximum normalized perveance per hole as long as similar nondimensional 
ized grid parameters are maintained. This result indicates that the 
screen hole plasma sheath characteristics would also scale with accel- 
erator system scale changes. 


APPARATUS AND PROCEDURE 
Sheath Probing Technique 

The screen hole plasma sheath was investigated at various grid 
geometry and accelerator system operating conditions. Figure 3 indi- 
cates the manner in which these data were obtained. In Fig, 3a a cross 
sectional view of the large seven hole two-grid accelerator system used 
for this study is depicted. Both the screen and accelerator grids were 
made from high purity sheet graphite. This material may be machined 
easily and is stable dimensionally at the high operating temperatures 
typical of this experiment 420°C). A seven hole hexagonal geometry 
was used for the accelerator system. With this aperture arrangement 
the central screen hole was surrounded by adjacent screen hole sheaths. 
In this manner, the central screen hole modelled a screen hole typical 
of a conventional ion accelerator system employing many thousands of 
apertures. The probing volume surrounding the central screen hole is 
shown by the dotted rectangle in Fig. 3a. The probe origin was set at 
the downstream face of the screen grid and the sheath was probed a 
distance of one screen hole diameter back into the discharge plasma 
from this point. The width of the probing region was 1.4 d^ (d^ = 
screen hole diameter) in order that sufficient data points could be 
taken to define the sheath overlap region adjacent to the screen grid 
webbing. 

Figure 3b indicates the manner in which the probe traversed the 
central screen hole. The sheath probe was set at each of the fifteen 
equidistant radial locations shown. At each location, the probe was 
swept axially through the sheath and the variation in local plasma 
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SHEATH PROBING 



b) 



Fig. 3 a) Screen hole sheath probing volume. 


b) Sheath probe location. 
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potential was recorded. In this way a full sheath profile, for a par- 
ticular accelerator system geometry and operating condition, was 
characterized by an array of local plasma potential values. 

Ion Source Design 

The accelerator system depicted in Fig. 3 was positioned over the 
masked down discharge chamber of an 8 cm diameter electron-bombardment 
ion source. This ion source had a mildly divergent magnetic field and 
used tungsten wire filaments as both the main and neutralizer cathode 
emitters. Argon propellant was used and all source operation was con- 
ducted in a 30-cm dia. pyrex bell jar. Average bell jar pressure was 
2.3 X 10-^^ torr. 

Ion source plasma conditions were determined from a small spherical 
Langmuir probe located inside the discharge chamber. This probe was 
offset slightly from the ion source axis and could be moved axially 
over a distance of some 3.0 cm within the discharge chamber. Plasma 
property measurements obtained with this probe were plasma potential 
Vp, electron density n^ and ion density n^. , Maxwellian electron tem- 
perature T|^ and the primary-to-Maxwell ian electron number density ratio 
Op/n^jj. Appendix A details the manner in which the spherical Langmuir 
probe traces were analysed to yield these parameters. 

Sheath Probe Design 

A floating emissive Langmuir probe^^ was used to study the screen 

I 

I hole plasma sheath and is shown in Fig. 4a. This probe had a simple 
dog leg shape so it could be positioned off the axis of the central 
screen hole thereby permitting the probe tip to move radially as the 
probe was rotated (Fig. 3b). An enlarged drawing of the probe tip is 
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Fig. 4 a) Sheath probe design details. 

b) Circuit schematic for floating potential 
measurements. 
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shown also to detail the basic construction of the probe filament, its 
pertinent dimensions and its attachment to the nickel support wires. 

In operation, the sheath probe is electrically isolated from the dis- 
charge chamber body and when the filament is heated to incandescence 
the filament floats at a potential very near the local plasma potential. 
This unique behavior of a floating emissive probe occurs for the 
following reason. Since the probe is floating, the net current at the 
filament surface must be zero. If the filament is cold (i.e., not 
electron emitting) this means an equal flux of plasma ions and elec- 
trons must arrive at the filament surface. However, because ions are 
so much heavier than electrons and their energy (for the electron bom- 
bardment discharge plasma under consideration) is less, they move con- 
siderable slower. Consequently, the initial electron flux to the 
floating filament surface greatly exceeds the ion flux. These electrons 
rapidly impress a negative potential upon the filament, impeding further 
electron flow until the equilibrium situation of equal ion and electron 
arrival rates is reached. The electron space charge so developed would 
ordinarily force the floating filament to assume a potential consider- 
ably negative of the surrounding plasma. The magnitude of this negative 
potential being of the order of the most energetic plasma electrons 
(tens of volts). If however, the filament is then heated to electron 
emitting temperatures the filament electrons so produced are accelerated 
away from the filament surface into the more positive plasma surround- 
ing it. In so doing, the filament electrons can neutralize the plasma 
electron space charge. This space charge neutralization is aided by 
increased filament electron emission with filament temperature increases. 
An equilbrium situation is quickly reached with the filament floating 
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very near plasma potential. Under these conditions, the flux of plasma 
electrons to the filament is countered by the outflux of filament elec- 
trons, which while possessing significantly lower energies can be pro- 
duced in much greater numbers (at least for the plasma densities under 
consideration in this experiment 'v. 10^^ m~^). Provided a sufficient 
filament electron temperature is maintained, to ensure adequate elec- 
tron emission, the emissive probe floating potential will follow 
closely local plasma potential changes. 

The instrumentation used to measure the probe floating potential 
is shown schematically in Fig. 4b. Briefly, an isolated battery supply 
provided heating power to the filament. The filament was connected 
through this battery supply to a high impedance electrometer set to 
measure voltage; the low side of the electrometer was referenced to 
screen grid potential. An emissive Langmuir probe characteristic show- 
ing the probe floating potential variation observed on the electrometer 
as a function of filament heater current is shown in Fig. 5. From 
Fig. 5 it can be seen that with the filament heater current between 1.1 
and 1.3 amperes the observed probe floating potential is not very dif- 
ferent from the actual plasma potential . Indeed, rather than take 
individual probe traces, the probe filament heater current was main- 
tained constant at 1.2 ampere and the local plasma potential variation 
recorded directly as the probe was moved through the screen hole plasma 
sheath. 

Sheath Probe Error 

Floating emissive Langmuir probes of the design shown in Fig. 4 
have an inherent error originating from the small D.C. voltage 


SHEATH PROBE FLOATING POTENTIAL-- volts 
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Fig. 5 Effect of filament heater current on the 
emissive sheath probe floating potential. 
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extending along the length of the heated filament. This potential drop 
was approximately 0.5 volt in this experiment. Another, more serious 
source of possible error when using constant battery supplied heating 
current, is the requirement of zero potential difference between the 
filament tip and electrometer input terminal. Figure 6 illustrates the 
circuit used to null this spurious voltage in this experiment. Briefly, 
the probe filament was heated by a 1.2 ampere current. The heated probe 
was then brought into contact with the discharge plasma (or upstream) 
side of the screen grid so that the incandescent filament was barely 
touching its surface. The nulling potentiometer depicted in Fig. 6 was 
then adjusted until the parallel resistance probe heating circuit was 
balanced and zero potential was observed between the filament tip and 
electrometer input. The sudden decrease in the electrometer potential 
reading when the heated filament tip touched the screen grid was also 
the manner by which a reference probe axial position was obtained. 

Both the nulling and probe position referencing technique were ex- 
tremely reliable and did not require continual adjustment. 

Another, more subtle, source of sheath probe error is associated 
with the magnitude of the double sheath which must surround a floating 
emissive probe. This double sheath forces an emissive probe to float 
somewhat negative of the surrounding plasma and is a consequence of 
the filament electron energy being considerably smaller than the plasma 
electron energy. Appendix B contains a theoretical derivation of the 
magnitude of this error, which suggests that a maximum error of one 
volt was to be expected. 

In practice, the screen hole sheath probe accuracy was checked by 
comparing plasma potential measurements obtained with it and with the 


SHEATH PROBE 
FILAMENT 



Fig. 6 Sheath probe voltage nulling circuit. 



19 


small spherical probe described earlier. The screen hole sheath probe 

consistently gave lower readings but the difference was less than S% of 

the spherical probe plasma potential readings in all cases. A more 

detailed discussion of the operating characteristics of floating emis- 

12 

sive Langmuir probes is presented by Kemp and Sellen. The position- 
ing error of the sheath probe was 0.004 d^ in the axial direction and 

t 0.02 d^ in the radial direction, 
s 

Sheath Contour Data Reduction 

The screen hole plasma sheath is really the set of equipotential 
lines that separate the region of homogeneous discharge plasma from the 
region of accelerated ions where there are no electrons. Equipotential 
contours describing this intermediate region were computer drawn from 
the array of data obtained after probing each screen hole sheath profile 
for each operating condition. First, each set of twenty axial sheath 
potential values, obtained directly from the sheath probe at each of 
the radial probing locations shown in Fig. 3b, was handled by a cubic 
spline interpolation routine. This computer routine produced fifty-one 
interpolated sheath potential values. Each of these points was 
separated by a distance of 0.02 d^ and together extended a distance of 
one screen hole diameter (1.0 d^) into the discharge chamber from the 
downstream screen grid origin. These fifty-one sheath potential values 
were curve fitted using a 10^^ order orthogonal polynomial regression 
analysis. Sheath potential values determined using the polynomial 
expression obtained from this regression analysis were then free of 
most of the scatter contained in the original data points. This 
scatter reduction was very important for electron and ion number density 
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calculations through the sheath. These calculations will be discussed 
in more detail later. 

Figure 7 is a plot of the twenty actual data points and the fifty- 
one points determined from the regression analysis polynomial for the 
sheath potential variation along the center line of a typical grid set 
for a standard operating condition. The agreement is excellent. The 
data in Fig. 7 has been normalized to a sheath potential one screen 
hole diameter back from the downstream screen grid face (the probing 
origin). The subsequent sheath potentials are negative to reflect the 
fact that the potential in the screen hole plasma sheath is less than 
that of the discharge plasma. Although interpolation routines were 
used for the axial sheath potential distribution, no such interpolation 
was done for the radial sheath potential variation when the sheath 
potential contours were plotted. 


SCREEN HOLE PLASMA SHEATH POTENTIAL- volts 
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RESULTS AND DISCUSSION 

Unless indicated otherwise, the following grid geometry and accel- 
erator system operating conditions may be assumed: 

center- to-center hole spacing, = 14.7 mm 

screen hole diameter, d^ = 12.7 mm 

total accelerating voltage, = 1100 volts 

net-to-total accelerating voltage ratio, R = 0.7 
screen grid thickness ratio, 

accelerator hole diameter ratio, d,/d = 0.64 

a s 

accelerator grid thickness ratio, t,/d = 0.37 . 

a S 

Effect of Plasma Density Variations 

Figure 8 shows the screen hole plasma sheath, represented as a set 
of equipotential contours, for the central screen hole of the large 
seven hole two-grid set shown in Fig. 3. The central screen hole in 
Fig. 8 is a cross sectional view and is drawn to scale as are the equi- 
potential contours. It should be noted that the path traced out by the 
sheath probe in Fig. 3b crosses the webbing of the adjacent screen 
holes at different locations. These locations were not always the same 
as those shown in Fig. 3b because each grid geometry tested was not 
orientated in exactly the same way about the center line of the ion 
source discharge chamber. Consequently, the path traced out by the 
probe, although always passing through the screen hole axis and being 
set at the same radial positions, did vary somewhat relative to the 
adjacent screen holes. To indicate this possible variation, the cross 
sectional view of the screen hole in Fig. 8 has the screen grid webbing 


NORMALIZED RADIAL POSITION ~ x/d 
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NP/H = 2.50x10-* amp/volt ®'2 
Bg/d, = 0.50 
Vd = 45.0 volts 

CONTOUR FROM 0.0 TO -25.0 volts 
CONTOUR INTERVAL OF 1.0 volt 



Fig. 8 Full screen hole sheath profile for standard 
grid set geometry and operating conditions. 
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extended out to the extremities of the probing volume, without regard 
to the actual location of the adjacent hole. 

The grid set used in Fig. 8 had an intermediate screen-to- 
accelerator grid separation ratio, of 0.50 and a discharge volt- 

age* Vjj, of 45.0 volts. This grid set was operated at a beam current 
or normalized perveance per hole value known to give the lowest ion 
beam divergence.® As mentioned previously, the sheath potential at a 
distance of one screen hole diameter back from the probing origin was 
defined as zero volts. Equipotential contours were plotted at one volt 
increments relative to this zero reference potential and are shown ex- 
tending to -25.0 volts (the plasma sheath potential being below that of 
the discharge, or bulk, plasma). The -25.0 volt contour line is very 
close to what was considered the sheath boundary. At this location the 
local electron density was only about 10% of the bulk plasma electron 
density. If the probe was moved much beyond this point no plasma 
electrons were detected and the probe accumulated a positive charge be- 
cause of the incident ion flux. 

The sheath surrounding the central screen hole of Fig. 8 is fairly 
symmetrical; the very slight asymmetry is believed to be due to the 
sheath probe entering different regions of the screen grid webbing at 
the extremities of its travel (Fig. 3b). The uncertainty in the axial 
location of the sheath potential contours was related directly to the 
spatial resolution which was 'V'0.02 d^ throughout these tests. Electric 
field vectors plotted normal to the equipotential contours give a fair 
estimation of actual ion trajectories (Appendix C). These vectors show 
how the average ion motion is directed through the screen holeasa conse- 
quence of the ion accelerating fields of the negative accelerator grid. 
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Perhaps the most important conclusion to draw from Fig. 8 is that the 
screen hole plasma sheath extends over a large distance, influencing 
ion and electron trajectories deep within the discharge chamber plasma. 

To illustrate this point, the Debye length for the plasma condition 
shown in Fig. 8 was 0.05 d^, whereas significant potential variations 
are evident over distances at least fifteen times this value. This 
result supports theoretical sheath width predictions of Grisham et. al.^ 
who have calculated potential distributions and ion trajectories for 
various ion accelerator systems. The work of Grisham et. al . and 
theoretical models developed by other authors will be discussed in 
greater detail in later sections of this report. 

Figure 9 makes a comparison between the sheath surrounding the 
central screen hole of the large seven hole two-grid set and the sheath 
surrounding one of the edge holes of this grid set. For clarity, only 
the -5.0, -15.0 and -25.0 volt contours are plotted in each case. Again, 
an intermediate grid separation was used and the grid set was operated 
at a normalized perveance per hole known to give the minimum beam di- 
vergence, Sheath distortion is evident for the edge hole. This dis- 
tortion tends to direct the initial ion trajectories to greater off axis 
angles, as evidenced by the electric field lines drawn in Fig. 9. Plasma 
density variations across the edge hole were primarily responsible for 
the sheath distortion shown in Fig. 9. An edge hole was probed also in 
a situation where the plasma density was known to be uniform. Under 

these circumstances, much less significant screen hole sheath distortion 

1 8 

was evident. It has been well documented that those ions emerging 
from the outer holes in a multiaperture accelerator system have very 
divergent trajectories and cause serious accelerator grid sputter 
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Fig. 9 Comparison of center and edge hole sheath potential 
contours (grid set edge towards vottom of figure). 
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erosion. The results of this study indicate that the decrease in plasma 

density as the discharge chamber edge is approached distorts the edge 

screen hole plasma sheaths giving rise to more divergent ion trajectories. 

The lack of adjacent screen hole sheath interactions for an edge hole 

does not significantly contribute to this sheath distortion and ion de- 

focusing effect. These results support earlier observations by 
4 

Kaufman. 

Effect of Beam Current 

13 6 7 

Previous experimental results ’ ’ ’ have shown that for a given 
total accelerating voltage ion beam focusing is poor at very low and 
very high beam currents. The best ion beam focusing occurs at a beam 
current intermediate of these two extremes. To investigate the mechan- 
ism of this focusing effect, the screen hole plasma sheath was probed 
during grid set operation at very low and very high normalized perveance 
per hole values. The screen hole sheath profiles determined from this 
probing are shown in Fig. 10 for the same grid set geometry as that used 
in Fig. 8. Here, the central screen hole of the large seven hole accel- 
erator system described previously was probed. Since the central screen 
hole sheath is fairly symmetrical (Fig. 8), only a half sheath profile 
is being presented. 

Figure 10 shows that at low normalized perveance per hole (or beam 
current) values the sheath boundary (i.e., the -25.0 volt contour line) 
is quite bowed and extends a considerable distance into the discharge 
chamber. By examining the electric field vectors, or ion trajectories, 
at the periphery of the sheath boundary, it is apparent that the ions 
emerging from this region have significant velocity components perpen- 
dicular to the screen hole axis. It is believed that these large 
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Fig. 10 Screen hole sheath movement and shape change with 
increasing beam current, or normalized perveance 
per hole, (half sheath profile). 
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off-axis ion velocity components produce the very divergent ion beam 
that is associated with grid set operation at low ion beam currents. 
Conversely, increasing the normalized perveance per hole, or beam cur- 
rent, by increasing the discharge chamber plasma density through in- 
creased cathode electron emission, moved the sheath closer to the 
screen hole and the sheath boundary became almost planar in shape. The 
trajectories of those ions emerging from the peripheral region of the 
sheath boundary (dashed arrow in Fig. 10) are now, if anything, directed 
slightly away from the screen hole axis; resulting in an ion beamlet 
that is fanning out at its extremities. It is believed that at this 
high beam current condition, such a large number of ions are being drawn 
through the accelerator system that electrostatic repulsion forces are 
fanning out the peripheral beamlet ions giving rise to poor ion beam 
focusing . 

The results of Fig. 10 indicate that the good focusing characteris- 
tics of beam currents intermediate of the two extremes shown are a result 
of a sheath boundary shape which gives minimum off-axis ion velocity 
components at an ion current density low enough to give negligible 
electrostatic repulsion effects. Figure 10 also shows that for the 
large normalized perveance per hole condition (dashed contour lines) 
the screen hole plasma sheath did not enter the screen hole. Previously, 
it had always been assumed that an ion accelerator system operating near 
its maximum normalized perveance per hole condition (as was the case for 
the dashed contour lines in Fig. 10) must have the screen hole plasma 
sheath located within the screen hole. A detailed experimental inves- 
tigation accounting for this behavior is presented in a later section. 
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Effect, of Grid Separation 

Figure 11 compares the effect on the screen hole plasma sheath as 
the separation distance between the screen and accelerator grids was 
varied. In all cases the grid sets were operated at a normalized 
perveance per hole value known to give the lowest ion beam divergence 
for that particular grid separation.® For clarity, only the -10.0 and 
-20.0 volt contour lines are plotted. Comparing the -20.0 volt contour 
lines shows that in all cases these lines were fairly coincident. Since 
the -20.0 volt contour is fairly close to the sheath boundary (Fig. 7) 
a tentative conclusion may be drawn. Namely, that to a fair approxima- 
tion the screen hole sheath boundary has the same position and shape for 
any screen-to-accelerator grid separation ratio when the grid set is 

operated at its minimum beam divergence condition. This result supports 

3 6 7 

earlier observations ’ ’ which indicated that the minimum ton beam 
divergence angle occurred at approximately the same normalized perveance 
per hole value irrespective of the screen-to-accelerator separation 
distance. Another feature of interest in Fig. 11 is that more plasma 
ions are directed away from the screen grid webbing and through the 
screen hole as the screen-to-accelerator grid separation ratio is re- 
duced. This effect is illustrated by the electric field vectors shown 
in Fig. 11. Deflection of plasma ions away from the screen grid webbing 
and through the screen hole as the grid separation is reduced would lead 
to increased maximum normalized perveance per hole values. Such an in- 
crease has been observed previously.® 

Effect of Discharge Chamber Parameters 

Tests were conducted to determine what effect discharge voltage, 
plasma electron temperature and primary-to-Maxwellian electron density 


NORMALIZED RADIAL POSITION*- x/d 


31 


NP/H = 2.45X id® amp/volt^ l?g/tls * 1.00 

NP/H = 2.50 X id® amp/ volt®^ C^/d, *0.50 

NP/H * 2.70X ld®amp/volt^2 Bg/d^ =0.25 


Vq* 45.0 volts 

CONTOUR FROM -10.0 TO -20.0volts 
CONTOUR INTERVAL OF 10.0 volts 



Fig. 11 Screen hole sheath movement and shape change with 

decreasing screen-to-accelerator grid separation, ^ , 

°s 

(half sheath profile). 
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ratio variations had on the screen hole plasma sheath characteristics. 
Figure 12 shows the results of these tests. Here, Fig. 12a depicts a 
half sheath profile for the standard grid set at typical operating con- 
ditions. Figures 12b and 12c show this same standard grid set operating 
at a higher discharge voltage and bell jar pressure respectively. In 
Figures 12b and 12c the sheath boundary was defined as that location 
where the local electron density was about 10% of the bulk plasma elec- 
tron density, (the determination of screen hole sheath electron density 
is presented In a later section). This definition was consistent with 
that used to define the sheath boundary in Fig. 12a. Comparing Figs. 12a 
and 12b it is evident that the sheath potential contours have similar 
shapes. However, the sheath boundary for a discharge voltage of 65.0 
volts is positioned farther from the screen grid than the standard 45.0 
volt discharge voltage condition of Fig. 12a. This displacement is 
non-uniform, with the sheath boundary adjacent to the screen grid 
webbing having been moved by a greater amount than the boundary position 

along the screen hole center line. Such a sheath movement with in- 

20 

creasing discharge voltage, had been theorized by Kaufman. The argu- 
ment being that adjacent to the screen grid webbing the plasma is shield- 
ing itself from a potential drop of the order of the discharge voltage. 
While, along the screen hole axis, the plasma is shielding itself from 
a potential drop of the order of the total accelerating voltage (Fig. 1) 
of which the discharge voltage is only a small fraction. Consequently, 
discharge voltage variations do not alter the total accelerating voltage 
significantly. Therefore, the plasma feels the accelerator system 
potential drop relatively unchanged and the sheath boundary is only 
slightly altered accordingly. 
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geometry and operating conditions (a), with sheath half 


profile for increased discharge voltage (b), and 
increased ion source operating pressure (c). 
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Of interest is to note that the sheath boundary shape associated 
with the high discharge voltage condition of Fig. 12b yields ion beam 
divergence angles lower by about 10% than those obtained at the standard 
operating condition of Fig. 12a. ^ It is believed that these lower ion 
divergence angles are a consequence of the sheath boundary being less 
bowed for the high discharge voltage than the boundary associa.ted with 
operation at a lower discharge voltage. This would mean that those ions 
emerging from the periphery of the sheath boundary in Fig. 12b would 
have smaller off-axis ion velocity components giving lower ion beam 
divergence angles than the ions emerging from the more bowed sheath 
boundary of Fig. 12a. Another interesting feature of Fig. 12b is that 
the sheath potential profiles adjacent to the screen grid webbing are 
quite flat and parallel the webbing surface. Consequently, no plasma 
ions are directed away from the screen grid webbing in this case. This 
is in contrast to operation at a lower discharge voltage, as shown in 
Fig. 12a. Here, the sheath potential contours adjacent to the screen 
grid webbing do direct some of the plasma ions away from the webbing 
surface and through the screen hole where they add to the extracted ion 
current. Figures 12a and 12b suggest that operation at a higher dis- 
charge voltage should result in lower maximum normalized perveance per 
hole values. This observation supports the result of earlier experi- 
ments® that showed higher normalized perveance per hole values were 
possible when the discharge voltage was decreased. 

Figure 12c shows the effect of increasing the bell jar background 
argon pressure by a factor of two over that used to obtain the results 
of Fig. 12a. From Fig. 12c it is evident that this pressure increase 
has reduced substantially the screen hole plasma sheath thickness. 


35 


With this large increase in neutral argon number density, the small 
spherical probe used to measure bulk plasma properties (Appendix A) in- 
dicated that the primary-to-Maxwell ian electron density ratio was approx- 
imately zero as compared to approximately 0.20 for the standard pressure 
condition of Fig. 12a, while the Maxwellian electron temperature was 
only slightly lower for the high pressure condition. It is believed 
that this near absence of primary electrons, with essentially only low 
energy Maxwellian electrons present, was responsible for the relative 
bunching of potential lines shown in Fig. 12c compared to those of 
Fig. 12a. Such a sheath thickness reduction would be expected to occur 
because the lower energy Maxwellian electrons would not penetrate far in- 
to the sheath. Consequently, the sheath potential gradients would be 
steeper. Figure 12c shows also that the sheath boundary for a high neutral 
argon number density has approximately the same position and shape as 
that associated with the standard operating condition (Fig. 12a). This 
similarity supports earl ier work that showed only slight ion beam divergence 
changes with argon gas flow, or bell jar pressure, variations.^ 

Effect of Screen Grid Thickness 

Screen grid thickness variations had a pronounced effect on the 
screen hole plasma sheath adjacent to the screen grid webbing. 

Figure 13 provides physical insight into these effects by showing how 
the screen hole plasma sheath moves as the screen grid thickness is 
varied. Over the central portion of the screen hole the plasma feels 
the negative accelerator grid potential most strongly; screen grid 
thickness and hence screen potential distribution changes have little 
effect in this region. Conversely, adjacent to the screen grid webbing 
the plasma sets up a potential distribution such that an equal ion and 
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electron flux arrive at the webbing surface. This latter effect is a 
surface phenomenon! only and is not dependent on the screen grid thickness. 
As the screen grid thickness is reduced, the plasma sheath boundary re- 
mains essentially fixed in position along the screen hole axis while 
adjacent to the webbing the sheath boundary follows the receding screen 
grid surface. However, with decreasing screen grid thickness the nega- 
tive accelerator grid potential is felt more strongly at the periphery 
of the screen hole. This tends to retard the downstream motion of the 
plasma ions and electrons after the receding screen grid in this region. 
The net result of these two competing effects is that eventually a 
limiting screen grid thickness is reached such that further thickness 
reductions do not alter appreciably the screen grid webbing sheath 
boundary position or shape. Although not shown in Fig. 13, other screen 
grid thicknesses were examined also. These tests indicated a limiting 
screen grid thickness ratio, t^/d^, of 'vO.05. 

Figure 13 shows that the motion of the sheath boundary after the 
receding screen grid surface gives rise to a sheath potential distri- 
bution at the screen grid webbing which tends to direct plasma ions away 
from the webbing and through the screen hole as screen grid thickness is 
reduced. This focusing effect is illustrated graphically by electric 
field vectors in Fig. 13. It should be noted that these ion trajectory 
vectors are similar to those presented in Fig. 11. Indeed, the web 
focusing effect evident in Fig. 11 was a consequence of physical pro- 
cesses similar to those described above. Except that in Fig. 11, de- 
creasing grid separation rather than decreasing screen grid thickness, 
resulted in the enhanced communication between the negative accelerator 
grid potential and discharge plasma at the screen hole perhiphery. 
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The arguments presented previously indicate that the position of 
the central portion of the screen hole plasma sheath would remain 
essentially unaltered with screen grid thickness changes. This be- 
havior is verified in Fig. 14 where screen hole sheath profiles are 
compared for a very thick screen grid " 0.49) and a conventional 

thin screen grid ("tg/dg = 0.18). Figure 14 shows that a plasma sheath 
has formed adjacent to the upstream surface of the thick screen grid 
webbing, but that the position of the central portion of the screen hole 
plasma sheath is still controlled by the negative accelerator grid po- 
tential and is relatively unchanged. The sheath probes shape. Fig. 4a, 
prevented the probe from entering very far into the screen hole and this 
is why only a portion of the -15.0 volt contour line and no -25.0 volt 
contour line is shown in Fig. 14, however, the trend is clear. Examina- 
tion of the sheath potential contours adjacent to the very thick screen 
grid webbing indicates that ions from within the screen hole could 
easily intercept the inner screen hole surface and recombine. Indeed, 
the normalized perveance per hole value indicated in Fig. 14, while 
only about 65% of the maximum value to be expected for the conventional 
thin screen grid (Fig. 2), corresponded to the maximum obtainable 
normalized perveance per hole for the very thick screen grid. 

The web focusing effect illustrated in Figs. 13 and 14 indicates 

that increased beam currents and decreased discharge chamber plasma 

losses (or increased beam ion production efficiencies) could be expected 

for screen grid thickness reductions. The trend of increased beam 

current, or normalized perveance per hole, has been observed previously 

o g 21 

for thin screen grids. ’ ’ 
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The different levels of ion recombination on the screen grid would 

★ 

be expected to be associated with different discharge loss levels. 
Figure 15 shows a plot of discharge loss, normalized to the discharge 
loss value for tg/dg ® 0.49, against screen grid thickness (the screen 
hole diameter d„ was held constant at 12.7 mm). This curve illustrates 
dramatically the adverse screen web focusing effect introduced with in- 
creasing screen grid thickness. 

The results presented in Figs. 13-15 are clear evidence that ion 

impingement on the screen grid webbing decreases significantly with 

screen grid thickness reductions. This has obvious implications for 

decreased screen grid sputter erosion and increased accelerator system 

22 

lifetime (a parameter of critical importance in ion thrusters ) as the 
screen grid thickness is reduced. In fact, it appears that an effort 
to thicken the screen grid to give longer erosion life may accelerate 
the erosion effect. Of course there are definite, but perhaps not in- 
surmountable problems in fabricating screen grids of the thicknesses 
suggested in Figs. 13-15. Also, without some sort of comparative 
erosion life test being conducted the lifetime conclusions presented 
here are tentative only. 

Effect of Screen Hole Shaping 

1 2 5 1 9 23 

Numerous workers ‘ * * have reported conflicting results as 

to the effect shaping the screen grid holes has on beam current and ion 
beam divergence. In order to further the understanding of the effects 


The parameter, discharge loss, is a measure of the beam ton production 
efficiency and gives a quantitative estimate of the energy (in eV) 
necessary to produce one beam ion. 
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of screen hole shape changes the screen hole plasma sheaths of variously 
shaped screen holes were examined experimentally. 

It was found that putting a 41® chamfer on the upstream (discharge 
plasma) screen hole face had only a very slight effect on the screen 
hole plasma sheath potential contours. Similarly, no change from the un- 
chamfered screen grid geometry was observed in the discharge loss level 

and maximum obtainable normalized perveance per hole. This null result 

23 

supports earlier evidence by Kerslake and Pawlik that screen grid hole 

chamfering is of little if any merit. 

1 5 

Some workers ’ have reported computer solutions and experimental 
results which indicate that counterboring the downstream screen hole 
face leads to significant screen hole sheath shape changes with lower 
divergence angles and increased beam current. Such a modification was 
made to the tapered screen grid discussed above. Figure 16 compares 
screen hole sheath profiles for this chamfered and counterbored screen 
grid to those of a conventional cylindrical screen hole geometry. Only 
slight differences in the sheath potential contours are observed. At 
the screen grid webbing the potential contours are fairly coincident. 
Towards the screen hole center more differences are apparent and the 
sheath boundary for the chamfered and counterbored screen hole (dashed 
-25.0 volt contour in Fig. 16) is displaced upstream slightly relative 
to the cylindrical screen hole sheath boundary. 

One would expect that by counterboring the screen hole the 
effective hole diameter has been enlarged. Such an enlargement must 
decrease slightly the positive potential that the screen grid webbing 
can impress at the screen hole center. Consequently, what shielding 
effect the screen grid could interpose between the discharge plasma and 
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Fig. 16 Effect of screen hole shape changes on sheath 
position and shape (half sheath profile). 
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negative accelerator grid will have been reduced. As a result, the 
plasma electrons would see a stronger retarding electric field as they 
approached the central portion of the screen hole. This stronger re- 
tarding field would prevent plasma electrons from penetrating through 
the ion acceleration region to a depth possible with the cylindrical 
hole geometry. The net effect would be to move the central portion of 
the screen hole sheath boundary slightly upstream, as is shown in 
Fig. 16. This trend agrees qualitatively with the previously mentioned 
computer solution predictions although it is much less pronounced than 
these solutions predict. 

The chamfered and counterbored screen grid shown in Fig. 16 gave a 
slight decrease in discharge loss (<3%) and a slight decrease in the 
maximum obtainable beam current {.< 3 %) compared to similar parameters 
recorded for the cylindrical hole geometry. No ion beam divergence data 
were obtained for any shaped screen hole geometry. Also, direct accel- 
erator grid impingement currents could not be measured accurately because 
of the relatively high facility background pressure ('v-2 x lO"** torr) . 
However, the slight screen hole plasma sheath shape and position changes 
apparent in Fig. 16 suggest beam divergence and direct ion impingement 
variations would be small as the screen hole shape was altered for thin 
screen grids. 

It is felt that the discrepancy between the results of the screen 
hole shaping tests presented here and those presented elsewhere^ *^*^*^^ 
lies in the screen grid thicknesses considered. As pointed out pre- 
viously, the screen hole plasma sheath is within the screen hole for 
large values of the screen grid thickness ratio only (Fig. 14). With 
the sheath inside the screen hole, it appears more likely that screen 
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hole shape changes would have a significant effect on the ensuing ion 
trajectories. References 1, 2, 5 and 19 all dealt with thick screen 
grids and large values of the screen grid thickness ratio In 

most cases these thick screen grids were a consequence of mechanical 
and thermal load requirements. From Figs. 13-15 it is evident that 
large screen grid thickness ratios must result in large ion source and 
screen grid thermal loads. It appears, then, that the degree of opti- 
mization for the accelerator systems studied in Refs. 1, 2, 5 and 19 
was significantly poorer than those typical of most ion thruster appli- 
cations; where thin screen grids are used. However, as evidenced from 
Figs. 13-15, even typical ion thruster accelerator systems (where 
^s^^s operating at higher than necessary discharge losses, 

if structural requirements are ignored. 

Experimental and Theoretical Comparison 

Figure 17 compares the screen hole plasma sheath obtained theo- 

4 

retically by Kaufman, using the space-charge-flow computer program of 

24 

Bogart and Richl.ey, against the sheath boundary (-25.0 volt contour 
line) obtained experimentally at the same grid geometry and operating 
conditions. Qualitatively, the sheath boundaries are similar but their 
positions are different. In the theoretical approach of Bogart and 
Richley the screen hole plasma sheath was defined as a single equipoten- 
tial surface at screen grid potential across which the electric field 
was set equal to zero. Consequently, the sheath obtained theoretically 
is represented as a discontinuity separating the discharge plasma from 
the region of accelerated ions, with the sheath terminating on the screen 
hole edge. This model was successful in predicting qualitatively many 
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Fig. 17 Comparison of experimental and theoretical 

sheath potential contours (half sheath profile). 
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4 

ion beam divergence effects. However, the present experimental inves- 
tigation indicates that the screen web focusing effects, screen hole 
sheath thickness and potential distributions are very important to a 
complete understanding of the ion extraction process. It is felt that 
the inability of the theoretical sheath model of Bogart and Richley to 
account for these effects is responsible for the poor quantitative 

4 

agreement betweeen theoretical ion beam divergence results and those 

obtained experimentally in previous studies 

Other, more rigorous, theoretical sheath models have been pursued 

but each has its own inherent difficulties and necessary approximations. 
25 

Wheal ton et. al . have listed the shortcomings of these various models 

while presenting the case for their own theoretical treatment, which is 

perhaps the most sophisticated model to date. The screen hole plasma 

sheath position and shape predictions obtained with this model agree 

qualitatively with some of the results of this study. However, the 

published results were for a large screen grid thickness ratio and it 

is not known how well the Wheal ton model could predict the screen web 

focusing effects found characteristic of thinner screen grids. Ion beam 

divergence predictions obtained with the Whealton model are contained 

in Refs. 5 and 25. These results contain a fair amount of scatter but 

do show trends qualitatively similar to those recorded by other 

workers.^ The accuracy with which the ion beam divergence 

data contained in Refs. 6 and 7 were obtained, coupled with the broad 

range of ion accelerator systems investigated and the success that has 
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been achieved by applying the results of this work, suggest that it 
should serve as the yardstick against which theoretical model predictions 
are measured. Similarly, the screen hole plasma sheath characteristics 
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presented and discussed in this section appear to define the important 
physical processes governing ion extraction from a discharge plasma. 

As such, these results should serve as a standard of comparison for 
those approximations of most validity in future screen hole plasma 
sheath models. 



SCREEN HOLE SHEATH CORRELATIONS 


Sheath Plasma Density Variation 

In the previous section, the screen hole plasma sheath was charac- 
terized by a set of equipotential contours. These contours defined the 
extent over which the plasma was perturbed by the accelerator system 
potentials and were very useful in estimating initial ion trajectories. 
Also useful is an understanding of the effect sheath potential gradients 
have on the local ion and electron densities. This information can be 
obtained directly from the sheath potential contours and the bulk plasma 
conditions. Briefly, a one-dimensional model was developed incorporating 
Poisson's equation, the ion energy equation and conservation of ion flux. 

The ions were assumed to enter the collisionless sheath region with the 
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modified Bohm velocity developed by Masek. * The second derivative 
of the potential with respect to distance was obtained numerically 
through the sheath from the experimental data. Using this as input to 
the model, it was possible to calculate the ion and electron density 
variation through the sheath. Appendix D details the development of 
this model . 

Figure 18 shows a plot of the ion and electron density variation 
through the centerline of the screen hole plasma sheath of the standard 
grid set at typical operating conditions (Fig. 8). The ion and elec- 
tron densities have been normalized with respect to their value a 
distance of one screen hole diameter back from the origin. The data 
presented in Fig. 18 were taken with two separate, but identically con- 
structed, sheath probes separated by a two month time period. What 
scatter is evident is believed lo be a function of the extreme 
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Fig. 18 Ion and electron density variation through a 


screen hole plasma sheath. 
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sensitivity in taking the second derivative of the sheath potential 
variation, rather than any inherent experimental inconsistancy. In 
fact, the sheath potential contours obtained with the different probes 
were virtually identical. Figure 18 shows clearly how the plasma ion 
density drops off slowly, as the ions are accelerated through the sheath, 
while in contrast, the plasma electron density drops off much more rap- 
idly due to the retarding field seen by the electrons. It should be 
noted that, to the author's knowledge. Fig. 18 contains the first ex- 
perimental plot of the ion and electron density variation through any 
plasma sheath. Some work was done by Goldan in the late 1960's on 
determining sheath potential variations using an electron beam probing 
method. However, he examined the sheath adjacent to a planar electrode 
and could only sense potential gradients an order of magnitude less than 
those detected during this experiment. 

Characteristically, all theoretical plasma sheath models assume a 
Maxwel 1 -Bol tzmann distribution of electron densities in the sheath region 
of the form 

"e " "o 

where V is the potential in the sheath and is negative, being set equal 
to zero in the bulk plasma, while n^ and T are the electron density 
and temperature (in eV) in the bulk plasma. In reality, the plasma 
electron population in most discharges is comprised of a Maxwellian 
and primary electron energy contribution.^ 

To test the validity of Eq. (1), the normalized electron density 
and sheath potential variation for the standard grid set geometry were 
plotted on a semi -log scale. Figure 19 shows the results of this 
plotting procedure. Here, the circled data points are for the standard 
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Fig. 19 Variation of electron density and potential 
through the screen hole sheath. 
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grid set geometry at typical operating conditions, the squares are for 
the same grid set geometry but twice the bell jar pressure (so that es- 
sentially only Maxwellian electrons were present), while the triangles 
are for grid set operation at a higher discharge voltage. These three 
operating conditions are the same as those discussed previously in 
Fig, 12. It is the normalized electron density versus sheath potential 
along the screen hole axis for the three sheath profiles of Fig. 12 that 
is being plotted in Fig. 19. From Fig. 19 it appears that the Boltzmann, 
or Barametric, equation (Eq. (1)) is a reasonable approximation to the 
observed functional form of the sheath electron density and potential 
variation. The effective electron temperatures determined from the 
slope of the lines of best fit, are consistent with the variations in 
the Maxwellian electron temperatures (Tni|), primary electron energies 
(<j>p) and the ratio of primary-to-Maxwell ian electron densities (Op/n^) 
of the different operating conditions. Calculations were made to 
determine if the effective electron temperature T, could be expressed 
quantitatively as a function of T^^, <|)p and n^/np. No successful corre- 
lations were obtained. It was felt that knowledge of the electron 
energy distribution function in the bulk plasma might yield an estimate 
of this effective temperature. However, this parameter is difficult to 
obtain experimentallyand was beyond the scope of this work. 

In summary, the results contained in Fig. 19 are experimental veri- 
fication of the validity of the Boltzmann, or Barametric, equation 

(Eq. (1)) for the description of electron densities in a plasma sheath. 

29 

Self, has commented on the range of applicability of Eq. (1) in a 
plasma sheath and these results support his theoretical predictions 
that any departure from Eq. (1) would be small. Although only the 
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electron density variation along the sheath axis has been presented in 
Fig. 19, similar results were also obtained for the electron density 
variation in the sheath adjacent to the screen grid webbing. 


Effective Screen Hole Sheath Area 

A method to determine an effective screen hole sheath area, as 
well as check the consistency of the entire experiment can be accom- 
plished by the following argument. 

Consider the ion current density equation given below 



"i 'i " 


( 2 ) 


where j Is the ion current density entering a screen hole, n^. and 

v^. are the ion density and ion velocity at any point through the screen 

hole sheath and e is the electronic charge. Setting as the 

effective screen hole sheath ion extraction area and making use of the 

1 6 27 

modified Bohm velocity developed by Masek ’ results in 




(3) 


Here, J is the ion current per screen hole, n^ is the plasma density 

at the screen hole sheath entrance, T„ and n /n„ are the Maxwellian 

m pm 

electron temperature and the primary to Maxwellian electron density 
ratio respectively, while K and M. are the Boltzmann constant and 
ion mass. 


Everything in Eq. (3) may be obtained experimentally except the 
effective screen hole sheath area A^, consequently, this parameter can 
be solved for directly. Numerous calculations were performed whereby 
the effective sheath area A^ was determined for a range of grid 
geometry and accelerator system operating conditions. The cumulative 
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experimental error associated with those measurable parameters in 
Eq. (3) prevented very accurate sheath area differences and trends to 
be obtained. However, for operating conditions where the screen hole 
plasma sheath had not entered the screen hole, application of Eq. (3) 
yielded effective screen hole sheath areas consistantly larger than the 
screen aperture area by not more than 20%. This result supports the 
observations of the previous section which clearly showed ion trajec- 
tories leaving a sheath boundary whose area was, in general, greater 
than the screen hole area. Perhaps the importance of the reasonable 
agreement between the calculated effective screen hole sheath area and 
that expected from the observed sheath potential contours is that the 
applicability of Eq. (3) has been substantiated. Inherent in the der- 
ivation of Eqs. (2) and (3) are the assumptions that the ion motion 
towards the screen hole plasma sheath is predominately one-dimensional 
with these ions assuming the modified Bohm velocity of Masek at the 
sheath entrance. These assumptions appear validated by the results 
discussed here. 


CONCLUDING REMARKS 


Results of the first comprehensive experimental investigation of 
the physical processes governing ion extraction from a plasma have been 
presented. The screen hole plasma sheath (.the transition region where- 
in significant ion acceleration and complete electron retardation 
occurs) has been defined successfully by equipotential plots for a 
variety of ion accelerator system geometries and operating conditions. 
These potential contours have shown that the screen hole plasma sheath 
extends over a large distance, influencing ion and electron trajectories 
significantly at least fifteen Debye lengths within the discharge 
chamber. The electron density variation within the screen hole plasma 
sheath satisfied a Maxwel 1 -Bol tzmann density distribution (i.e., the 
Barometric equation) at an effective electron temperature dependent on 
the discharge plasma primary to Maxwellian electron density ratio. 
Similarly, plasma ion flow up to and through the sheath was predominately 
one-dimensional and the ions entered the sheath region with a modified 
Bohm velocity. Only at large values of the screen grid thickness ratio 
did the screen hole plasma sheath enter the screen hole. However, a 
significant screen webbing ion focusing effect indicated poor plasma 
generation efficiency and low extracted ion currents with large values 
of the screen grid thickness ratio. 

The results of this work define those parameters of most importance 
to the ion extraction process from a plasma. It is hoped that the 
development of successful theoretical Ion extraction and focusing 
models might be aided by the physical insight embodied within this work. 


APPENDIX A 


THICK SHEATH SPHERICAL PROBE ANALYSIS 

The small spherical Langmuir probe used to measure discharge plasma 
properties is shown in Fig. A-la. Plasma conditions in the ion source 
discharge chamber used for the screen hole sheath study were such that 
the Debye length (the charged particle shielding distance) at typical 
operating conditions was 0.5 mm or greater. This large Debye length 
necessitated the use of a thick sheath probe analysis. A typical 
spherical probe trace is shown in Fig. A-lb. This probe trace was 
analyzed in the following way. 

In the retarding region the electron current to the probe was made 
up of thermalized electrons in a Maxwellian energy distribution and 
unthermal ized primary electrons of energy near the plasma discharge 

13 

voltage. This electron energy distribution has been studied by Martin 
for an argon discharge. From Beattie^^ the electron current to the 
probe in the retarding region of the probe trace shown in Fig. A-lb is 
given by 

I = Bi + B2 V + B3 exp(B4 V). (A- 1 ) 

Here, the primary electron current is the linear portion while the ex- 
ponential portion pertains to the Maxwellian electron current. A 
non-linear least squares fit technique was used to fit Eq. (A-1) to 
the retarding portion of the probe trace depicted in Fig. A-lb. This 
procedure was a computer library routine and converged rapidly to give 
the coefficients Bj, B 2 , B 3 and B 4 . Wrth Bj and B 2 known the straight 
line primary electron current contribution was subtracted off the entire 


S8 




b) 

Fig. A-1 (a) Spherical Langmuir probe design details. 

(b) Typical thick sheath spherical Langmuir probe trace. 
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Langmutr probe characteristic shown in Fig. A-lb. The equation to the 
electron retarding region now described Maxwellian electrons only and 
became 


Ln I = LnBa + B 4 V. 


(A-2) 


Also, the Maxwellian electron current in the saturation region of the 

1 5 

probe trace is, for a thick sheath spherical probe, 

I = B 5 V + Be . (A-3) 


Where B 5 is the slope of the electron saturation region shown in 
Fig. A-lb after the linear primary electron current contribution has 
been subtracted off the probe trace. The magnitude of this slope is 
related to the conductance of the plasma. Similarly, 65 may be obtained 
by solving Eq. (A-3) at any current-voltage point in the Maxwellian 
electron saturation region. Substituting Eq. (A-3) into Eq. (A-2) gives 

B5 B5 

I = (Be - ^ Ln B 3 ) + ^ Ln I . (A-4) 


Equation (A-4) was solved iteratively to obtain the Maxwellian electron 

saturation current 1^^^ and then either Eq. (A-2) or (A-3) was used to 

evaluate the plasma potential V . Where I . and V are the current 

p sat p 

and voltage values at that critical point in the Langmuir probe trace 

where plasma electrons are neither retarded or attracted by the probe. 

Using these results plus the inverse of coefficient B 4 (which is the 

Maxwellian electron temperature T^^ in eV) the electron density n^ was 

1 5 

detennined using the following equations. In the retarding region: 


n = 2.969 X lO^^ 


sat 


m 


m' 


_3 


(A-5) 
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In the saturation region: 

/r 

n^ = 2.969 X . (A-6) 

e r2 

P 

Here, rp is the spherical probe radius in meters. The electron density 
values computed from Eqs, (A-5) and (A-6) agreed within 20% consistently. 
The average of both values was used to determine discharge plasma elec- 
tron and ion densities. Also, the ratio of primary to Maxwellian 

electron number density n /n^^ was evaluated for each probe trace. The 

14 

equation used in this determination is given below 

!e . 2_ ( _!"! )■" ^ . (A-7) 

"m Vp + B 1 /B 2 sat 

Where I . is the primary electron current reaching the probe at 
prim 

plasma potential and was calculated by using the linear portion of 
Eq. (A-1). While the quantity, Vp + B 1 /B 2 , is equal to the primary 
electron energy, <j>p in eV. 



APPENDIX B 


SHEATH PROBE ERROR 

An emissive Langmuir probe has been used in the screen hole plasma 
sheath study. When surrounded by a plasma, and with the filament heated 
to incandescence, an emissive probe will float at a potential near the 
local potential of the surrounding plasma. The error, or the difference 
between the emissive probe floating potential and plasma potential is the 
result of two effects. The first is the potential variation along the 
hot filament wire because of the ohmic heating voltage drop. For this 
experiment a carefully balanced D.C. battery supply was used to heat the 
filament. The voltage drop along the length of exposed heated filament 
has been estimated to be the order of 0.5 volt. The second source of 
error is due to the formation of a double sheath around the floating 
probe. This double sheath Is a direct consequence of having more elec- 
trons produced by the hot filament than are needed to satisfy the re- 
quirement of zero net current between the floating probe and plasma. A 
quantitative estimate of the voltage drop in this double sheath is the 
purpose of this appendix. 

For this analysis the simplified case of an infinite planar elec- 
tron emitter is considered and the analysis is one-dimensional. The 
filament is floating and assumed to emit more electrons than necessary 
to satisfy the zero net current requirement (this represents a condition 
of operation typical of the emissive probe used in the screen hole 
plasma sheath study). The double sheath thus formed is shown as 
follows. 
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Define: 


'M 


Plasma potential . 

Planar filament floating potential. 

Double sheath potential minimum. 

Potential at any point between filament and plasma. 


Po-j = Bulk plasma ion charge density. 

Pog = Bulk plasma electron charge density, 

Pof = Filament electron charge density at filament surface. 

p-» P » Pf = Plasma ion, electron and filament electron charge density 

’ ® at location where potential is V, 

Pm» p-«, Pjtm “ Plasma ion, electron and filament electron charge density 

n em rm V = V,^,. 


To estimate the potential difference it is necessary to solve 
Poissons’ equation, which in one-dimensional form is; 
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Multiplying both sides by dV/dx and integrating with respect to X gives 



At the double sheath boundaries the voltage gradient is zero. This 
leads to the boundary conditions: 


dx 


= 0 — 


V=V^ 


dv 

dx 


= 0 . 


v=v_ 


0 m 

Applying these boundary conditions to the equation above yields 



0 . 


(B-1 ) 


Also, since the planar filament is floating, the net current between 
the filament and bulk plasma must be zero at any point. Examining the 
currents that must be present at the double sheath potential minimum 
results in the following equation: 


"fM ''fM * "IH “,-H ' "eH 

where and are the velocities of the filament electrons and 

plasma ions and electrons at the potential minimum. 

Equations (B-1) and (B-2) need now to be solved simultaneously to 
obtain the potential difference First, Eq. (B-1) is constructed 

piecemeal by considering each of the species present. 


Plasma Ions 

The plasma ions are assumed to enter the double sheath region from 
the plasma with the Bohm velocity,^® since no primary electrons are 
being considered. These ions are accelerated as they travel from Vq to 
and are retarded slightly as they go from Vj^ to V^. Restricting 
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ourselves to the accelerating region only, the Ion energy equation 
can be written as follows: 


2 kTj^ - 2 V? - eCV^ - V) 


where is the temperature of the Maxwellian plasma electrons and v^. 
is the ion velocity at any point between Vq and V^. Re-arranging gives 


"i = 


kT, 


H * (V. . V) 


From conservation of ion current a second equation is obtained 


kT, 


M 


h 


*^oi \ m. j ” ^i '^i 


Substituting the above form of v. gives: 

Pi Poi ~ 

Integrating this expression with respect to V gives 
V V 

/ n ^ lf„ (''o 


m 


m 


setting 


9 

t2 = 1 + ^ (Vq - V) and dV = - t dt 
M 


and noting that when 


V = V„ , t = 1 


and when 


one obtains 




V = V„ , t = [1 + If- (V, - v„)] 


M 

kT. 


■M 


r'^o kT„ 

I Pi'iv = --rp„i 

\ [1 


dt 

2e 
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or Vq 

/ '> 1 '"' " ^ I ^ I • ( 8 - 3 ) 

Vh 

Plasma Electrons 

Electrons from the plasma see a potential hill as they enter the 
region between and and are retarded. This potential hill, as 
viewed by the plasma electrons, is shown below. 



Consider some point x' lying between and The plasma electron 
space charge at x' is due to two electron velocity groups. Plasma 
electrons with initial velocities between (V^-V)]* and (V -Vn.)]^ 
reach and go past x', are reflected back from the potential hill and 
pass through x' again on their way back into the plasma. Hence, these 
electrons contribute twice to the electron space charge at x'.^^ Also 
those electrons with initial velocity between ~ 
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through x' and contribute to the charge density at this point. But 
these electrons cross over the potential hill and never return to 
the plasma The plasma electron charge density at x' then takes the 
following form. 

r|e /y _y \-h 

e 

Oe = 2^oe/ * “-oe / 

e e 
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Similarly, substituting energy equation 

e 

gives: 


"(''o-'m) ' i ”e''e * 


or 


''I “ • 


Hence, the expression for becomes: 

- Cr(v-v„)]’^ 

■e(V^-V) ' ' '*® 


/ ^ 

Pe =Poe (,2^ ) 


kT, 


M 


/■ -.-i 


2 |exp 2 kj^ dVg + I exp 


a 

f 


■My 


2kT, 


av 


M 




One can now set t 


and note that v = 
e 


I^Y 

\2kT„ J 

"ft 


'2 / \h 


"e =[m, 


dv. 


M 


corresponds to t = -r|- (V-V„) 

. M 


Substituting these results back into yields: 


^e = Poe 


(2,kT„)“PUT„ A.H^ ; 


2jexp{-t^)6t +fexp(-t^ 
0 re 




t^)dt 

h 


But, yexp(-t^)dt = erf(x) 




exp(-t )dt = f [l-erf(x)] 


and 


CM <U 
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therefore 

Pe " I Poe ®^P (lT;;j^ ) * 


Integrating this expression with respect to V gives 



Unfortunately this integral is too cumbersome to be evaluated with ease 
analytically, if at all. However, it may be approximated fairly 
accurately by the use of Simpsons' three point rule. The result is 
given below: 



M 


+ 4 exp 



(B-4) 
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Filament Electrons 

Electrons from the filament see a potential hill as they leave 
the filament. 



Only those filament electrons with an initial velocity between 

[ r «*> can contribute to the filament electron charge 

e 

density at x'. This charge density is given by the following expression 


oo 


P4: = p 


fH 

[r <W]^ 


Where f( is the one-dimensional Maxwellian velocity distribution of 
electrons at the filament surface 





-%Jk] 

2kTf ) 
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and is the filament electron temperature. Using the filament 
electron energy equation 


I “ ? Vf + ® 


the Maxwellian velocity distribution function at any velocity v^ is; 


ml 


""pl- 2 Frr 


-e(V^ - V) 


When v^^=^ (V^ - the preceeding energy equation gives 

e 


e(Vf - V =|Mgv2 + e(V^ - V) 

''f ‘ [c - ''h>] • 


"e 


Substituting these results into the expression for gives 

/ «e \>‘‘ - ''V” ... 

'’f '’fM \^2itkTj, j \ kT^ )J \ 2kT 


/-VfV 

\ 2kT^/^''f 


If one sets t = 


K /2kJ.\h 

2 \a^j Vf and dv^=^-pj--j dt 


and notes that v^ = 
corresponds to t = 


one obtains = PfMl'^ 




-e(V,-V)\/2kT.\!5/'” 
kl! )(-H- J exp(-t2)dt 


[kTf*''"''H’] ^ 


or 
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Pf 2 PfM 


-e(Vf-V)' 

~TjZ 


1 - erf I 


e(V-V„)U 

kT, y 


Integrating both sides with respect to V gives 


1 

/ 


'M 


"o dv - 1 

PfdV - 2 


/e(V-VJ 


1 - erf 


'e(V-V^) 


'M 


Typically, |v,^-Vf| ~ ~ 0.2 

and jvj ~ one volt or more. 


eV 


kT. 


dV. 


Hence, a good approximation is to set 
e(V-V„) 

If one sets t^ = ” 


kT. 


then 2tdt = -j^ dV 


When V = Vjjj, t = 0, also when V = Vq, t = j^-j^ ^ * 

Substituting these results back into the previous integral gives 
M \i 


e(V-V„) 


. f o , Ao _ /e(V- 

J Pf")'' ’■zomj «P 


dV 


1 

- o P 


2 "^fM e 


[kT^ '^Om] 

-J" exp(t2) erf(t)dV 


0 




Now for t^ > 5, in this instance a very good approximation, one can use 
the following simplification 


exp(t^) erf(t) = — 


1 


1 


If* (t + 2^) 


Thus the above integral becomes: 
V 


r 


1 


PfdV = p 


fM 


& [». {^) ■' 1 J 


'M 
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2kT, ^ rIt ''oh 


- P 


fM 


r Kl^ UM ^ 

/n e J ^ 


dt 


But 

— n — dt = / -T-l- — „ dt = t — - arctan 

t+Jt ^ 


(t/n 


and therefore 
V 


/ 


'M 


PfdV 


^fM 


r kT r r ® 1 r -I ^ 

[21- [exp(-i?r; ’ "'j ■ ^ [107 ''oh] 


kT 

+ i arctan 

e(2ir)'^ 


r^v 

L kT^ ''OM J J • 


(B-5) 


Substituting Eqs. (B-3), (B-4) and (B-5) into (B-1) and remembering 
that Pg and are negative, gives: 



Now Eq. (B-2) had the following form 


•^fM ''fM ■*■ ^iM ''iM " ^eM'^eM * 

From the previously calculated charge density expressions and initial 
velocity assumptions each term in this equation can now be determined. 
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^fM PfM 


fM 


v]' 

% ''oh ] 


^•M " ^01 


=r^+^v 

1M Lm^. M. ''omJ 


_ 1 '^OM 

^eM - 2 %e 


kT, 


M 


eM 


= [!l!!tl.2e^ ]% 
[ttM^ ''OM I 
'- e e -* 


Substituting back into Eq. (B-2) and dividing by gives; 

kT. 


8kT 

t:M 


f . ^ 


P 


01 


'M . 2e 


M^. M^. ''OM 


1 ^oe 
= 2 

PfM 


*M *f'J 

PfM 

1 

L kT 

t" ''oH'i 


2e y 


ttM 

L e 

M^ ''OM 
e 


M 


or 


i!oe 
oi _L 2 PfM 


-e V 


exp 


OM 


kT, 


'fM 


M 


8kT„ . 

— \i 

irM^ '^OM 

e e 


— t. + ( V -V i 

irM^ M^ ^''m ''r 

e e 


M^. M^ ''OM 

' * lf„ ''oM, 


Dividing both sides of Eq. (B-6) by and substituting in the above 
p. 


expression for yields, after re-arrangement, the following equation; 
PfM 
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Standard operating conditions of the ion source used in the screen hole 
plasma sheath study resulted in a Maxwellian electron temperature of 
= 8 eV (92752°K). While the probe filament electron temperature was 
approximately = 0.2 eV (231 9°K). Substituting these results in the 
above equation gives: 


1 + 


0.25 0,5 ^exp(-0.125 .582 x1012_3.516x10HVqj^]^ 


- 1 .386x105 


1 .91 6x1 0^+4.790x1 O^VoM 


- 


1 + 0.25 V 


OM 


1 ^oe 


12 p 


fM 


[exp(-0.125 Vq^) +1 + erf (0.1 25 ) + 4 exp( -0.0625 Vq^) 

[1 + erf (0.0625 

= 0.1[6Xp(5Vq^)-1] -0.113[5Vom]'^ + 0.080 arctan[10V^M]^ . (B-8) 


'OM 


OM- 


0 G 

Solutions to Eq. (B-7) were sought for different values. The 

PfM 

results are tabulated below. 
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Table B-1 Double Sheath Potential 


‘^oe 

'’fM 

Vqm - volts 

0.18 

0.6 

1 .0 

1 .5 

10 

2.0 

100 

2.5 

1000 

3.0 

10000 

3.45 

100000 

3.92 



00 

No solutions were possible for values below 0.18. The formalism 

PfM 

of Eq. (B-7) seems to suggest that this lack of solution was because a 
greater plasma electron current was being requested than the initial 
conditions of the problem permitted. With decreasing filament electron 
emission (i.e., increasing ^ — values) Vgj^ increases slowly. Since 
the model has assumed that there are no primary electrons present, the 
cold filament would be expected to float several volts below plasma 
potential. This result is born out by the trend in the above table. 

The results of the model indicate that if the sheath probe filament 

Pop 

were emitting electrons sufficient to make the ratio — 1 then the 
probe potential measurements in the bulk plasma would be in error by 
about one volt. Since the plasma electron temperature used in Eq. (B-8) 
corresponded to a plasma potential of about 43 volt this means the 
screen hole sheath probe error is about 2-3% of the true plasma potential 
and on the low side. 


/ 
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Filament Electron Emission 

The preceeding model has assumed that more electrons are being 
produced by the filament than are necessary to maintain zero net current 
between the filament and bulk plasma. A worse case is to assume that 
all the filament electrons produced are used to maintain this zero net 
current requirement. Under these conditions the electrons leave the 
filament with space charge limited emission. Consequently, no virtual 
cathode is formed and only a single sheath is present. Hence, may 
be replaced by the filament electron charge density at the filament. 
The filament temperature, under typical operating conditions, has been 
estimated to be 2750°K. This corresponds to a filament electron tem- 
perature of 0.24 eV. Using this temperature and the Richardson-Dushman 

equation, the charge density of 8 eV plasma electrons to the charge 

^oe 

density of filament electrons, , is about 0.1. This result indi- 

Pfo 

cates that the emissive probe filament is more than capable of providing 
the emission levels to keep the magnitude of the double sheath potential 
hill to acceptably small values. 

A flat planar filament has been assumed for the sheath probe error 
model. In reality, the filament was a thin hairpin of tungsten wire. 
(Fig. 4a). However, at typical discharge plasma conditions the thick- 
ness of the double sheath surrounding the probe filament has been 
estimated to be at least of the order of the filament wire diameter. 
Under such conditions, a flat planar filament is a fair approximation 
to the actual physical situation. 

Probe Error in Screen Hole Sheath 

When the emissive filament enters the screen hole plasma sheath 
the double sheath adjacent to the probe changes. The ions are being 
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accelerated through the screen hole sheath but the ion current to the 
filament remains constant from the ion flux conservation requirement. 
Also, the filament heating power and hence temperature remains fixed 
and consequently the number and energy of filament electrons available 
for electron emission is unaltered. What does change however, is the 
filament directed plasma electron current at any point through the 
sheath. This current decreases with distance through the sheath while 
the temperature of this plasma electron current remains unchanged. 
Therefore, one would expect a decreasing filament emission current the 
farther in the sheath is probed in order to maintain zero net current 
between the quasi plasma in the sheath and the filament. If the probe 
were pushed through the sheath so that no significant number of plasma 
electrons were present, the probe filament would be expected to accumu- 
late a positive charge due to the incident ion flux. This behavior 
was observed experimentally. 

The model developed here has taken no account of the presence of 
primary electrons in the discharge plasma. This assumption would seem 
valid because the ratio of primary to Maxwellian plasma electrons is 
approximately 20% for the discharge plasma under consideration. The 
importance of primary electrons, as far as the probe error is concerned, 
is that they are the specie which defines the downstream screen hole 
sheath boundary. Since the primary electron number density is rela- 
tively small those primary electrons which define the sheath boundary 
contribute a small electron current to the probe which drops off 
rapidly as the sheath boundary is approached. The definite size of 
the sheath probe filament precludes adequate resolution of this termi- 
nation region. This resolution loss becomes significant only when the 
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sheath has been probed to a depth where the sheath potential has dropped 
to a value about 40% of the plasma potential. At this point, the po- 
tential is decreasing with distance so quickly that two filament widths 
account for the remaining 40% of plasma potential (Fig. 7). 


APPENDIX C 


PLASMA ION SHEATH TRAJECTORIES 

It IS of interest to know how faithfully plasma ion trajectories 
follow the electric field lines in the screen hole plasma sheath. 
These ion trajectories may be determined by considering the following 
diagram which shows a portion of a typical screen hole plasma sheath. 


^2 

Ar 

AZ 


Here, Vj and V2 are the magnitudes of any two equipotential contours, 
E(aS) is the average electric field strength between these contour 
lines and is a function of their separation AS, while E^(Ar) and 
E^CaZ) are the radial and axial components respectively of E(aS). The 
acceleration a plasma ion receives upon passing between Vj and V2 has 
the following components 
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E (Ar)e 
”i 

E7(AZ)e 

3z 

"i 


where and e are the ion mass and electronic charge respectively. 
Equations of motion describing ion motion between Vj and V 2 are given 


+ 2a Ar 
r or r 


(C-3) 


v^=vJz+2a^AZ . 


(C-4) 


Where v^ and v^ are the radial and axial ion velocity components 
at V 2 while v^^ and v^^ are the initial radial and axial velocity 
components at Vi. Substituting Eqs. (C-1) and (C-2) into Eqs. (C-3) 
and (C-4) respectively gives 


2eAr E„(Ar) 

v2 = v2 + — -r 

r or M. 

2eAZ E,(az) 

v2 = v2 + ? 

z oz 


(C-5) 


(C-6) 


Equations (C-5) and (C-6) were used to calculate plasma ion trajectories 
through a typical screen hole plasma sheath profile. Ions were assumed 
to enter the screen- hole plasma sheath at a distance of one screen hole 
diameter from the origin with the modified Bohm velocity developed by 
Masek.^^*^^ This velocity is given by 


V 


(C-7) 
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where e is the electronic charge and T and are the Maxwellian 

m 

electron temperature and primary-to-Maxwellian electron density ratio; 
these parameters were equal to 7.8 eV and 0.20 respectively for the 
screen hole sheath used in the ion trajectory calculations. Ion 
trajectories were investigated for initial ion velocities parallel 
to and 10° above and below the axial co-ordinate direction. In the 
case of the non-axial initial velocities, the appropriate cosine and 
sine were multiplied by Eq. (C-7) to define the initial axial and radial 
velocity components respectively. By calculating the angle between the 
axial and radial velocity components obtained after each Vi V 2 poten- 
tial step, the overall ion trajectory through the sheath could be 
plotted. 

Figure (C-1) shows the results of these ion trajectory calculations. 
Ions entering the sheath parallel to the sheath electric field vector 
lag behind this vector, but not significantly. Similarly, those ions 
entering the sheath with the off-axis angles shown are eventually brought 
into reasonable alignment with the local sheath electric field vector. 

It should be mentioned that plasma ions assume the modified Bohm velo- 
city as a result of the slight potential gradients which extend back 
into the center of the ion production region of the discharge chamber. 
Since plasma ions are created with a randomly directed energy equal to 
the discharge chamber temperature ('vO.06 eV), they must become aligned 
rapidly with the bulk plasma electric fields, which although small, can 
and do accelerate these ions up to a velocity of several eV before they 
enter the screen hole plasma sheath. Consequently, one would expect 
most of the ions to enter the sheath parallel to the sheath electric 
field vector. For an ion to have even a 10° departure from axial 
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ELECTRIC FIELD VECTOR 

CALCULATED ION TRAJECTORIES 

CASE A : ION ENTERING SHEATH PARALLEL 
TO ELECTRIC FIELD VECTOR 

CASE B : ION ENTERING SHEATH 10® ABOVE 
ELECTRIC FIELD VECTOR 

CASE C : ION ENTERING SHEATH 10® BELOW 
ELECTRIC FIELD VECTOR 



Fig. C-1 Ion trajectories through a screen hole plasma sheath. 
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alignment would be unlikely because the ion production rate near the 

22 

screen grid for the ion source used is small. Also, the probability 
of ion-neutral atom collisions or charge exchange processes was small 
because of the low ion source propellant pressures used. 

In summary, it may be stated that ion trajectories through the 
screen hole plasma sheath follow the local sheath electric field vectors 
to a fair approximation. 


APPENDIX D 


SHEATH ION AND ELECTRON DENSITY VARIATION 


In order to detemine the electron and ion number density varia- 
tion through the screen hole plasma sheath the three equations follow- 
ing were considered: 


v2y = 


e(ng-n^) 


(D-1) 


n.v. 
1 1 



(D-2) 


jH.Vf 


KT 


n 


ri 2 ■■■i'- M.“ ^ • 

1 m 


(D-3) 


Here, Eq. (D-1) is Poisson's equation for electrons and ions where V 

is the sheath potential and is negative. Equation (D-2) equates the 

ion flux at any point in the sheath to that at the start of the sheath 

1 fi 71 

where the modified Bohm velocity criterion of Masek ’ is assumed to 

apply. Equation (D-3) equates the ion energy at any point in the 

sheath to the ion kinetic energy at the sheath entrance plus that 

gained as the ions are accelerated through the sheath. As before, n^ 

and n. are electron and ion number densities at any point in the sheath, 

n 

T|^ is the Maxwellian electron temperature, ^ is the primary- to- 

"m 

Maxwellian electron density ratio and M., e, K and are the ion mass, 
electron charge, Boltzman constant and free space permitivity 
respectively. 

Azimuthal symmetry is assumed and consequently a one-dimensional 
analysis is pursued. In this instance the co-ordinate direction of 
interest is perpendicular to the screen hole sheath centerline. With 
these assumptions Eq. (D-1) reduces to 
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e(n^-n,) 


Rearranging Eq. (D-3) gives 

KT n„ 1 

V = (r — ni (1 + _E )i _ ±®V 

i M. n„ M. ^ • 

1 m 1 

Substituting this expression back into Eq. (D-2) yields 


nT-jP- (1 + -^ )]' 

0 M . n„ 

] m 

KT n o If 1 

(r—El (1 + _E )] _ )'^ 

M. n„ M. ^ 

1 m 1 


This expression may be substituted back into Eq. (D-4) to obtain a 


relation for n that is 
e 


KT„ n , 

-oC IT ^ >3^ 

] m 

£([-/(! - ^) 


KT n 1 

n [ ncr (1 

[o . ° "m 

e az^ KT^ n„ i 

(f —JO. (1 + _B. )•] _ 

M,- M, ' 

1 m T 


Equations (D-5) and (D-6) may be solved numerically to find n. and 
n^, since all other parameters can be obtained experimentally. 
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